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We consider two component Bose-Einstein condensates prepared in a single spin projection. In 
that case, density and spin waves can be described as distinct elementary excitations. This gives 
rise to separate critical velocities for the breakdown of the motional and the spin superfluidity. We 
demonstrate this behavior considering the interaction of the condensate with both a mechanical and 
a magnetic defect and calculate the corresponding drag forces. 
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Introduction. Superfluidity [lj is undoubtedly the 
most striking and famous phenomenon linked with the 
Bosc-Einstcin condensation (BEC), as it corresponds to 
regimes in which the condensate is stable against the cre- 
ation of excitations. According to the Landau criterion, 
it occurs whenever the condensate propagates at a ve- 
locity smaller than the speed of its excitations The 
latter are density waves, or bogolons 0], characterized by 
a linear dispersion for long wavelengths and a density- 
dependent celerity of sound. When the internal spin de- 
gree of freedom of the particles is taken into account, 
the situation becomes even more interesting. Indeed, 
spmor BECs 0, [|, which are not restricted to atomic 
systems but also include other systems such as exciton- 
polaritons @, Q and magnons 0, Q , have already brought 
remarkable discoveries, such as the formation of topolog- 
ical defects, exotic phases and spin textures (lollllj. The 
physics of spinor BECs becomes even more reach in the 
presence of spin-orbit interactions fl^4l6j . 

In general, the superfluidity criteria for spinor BECs 
are less trivial, since the critical velocity is not uniquely 
defined 0] . As a consequence, spin and density degrees of 
freedom are mixed and the channel associated with pure 
spin excitations [l7l - fl9j , in analogy to magnetism ptj ] , is 
harder to identify. Therefore, the task of isolating, pro- 
tecting, and controlling pure spin excitations can be of 
great interest for applications. The study of such sys- 
tems leads to the understanding of spin superfluidity as 
a phenomenon distinct from the usual density (Landau) 
superfluidity: the onset of a spin-polarized (or magne- 
tized) flow without spin wave excitation. In fact, spinor 
BECs possess undeniable similarities with magnetic sys- 
tems: depending on the relative strength between the 
intra and inter-component interactions (tunable e.g. via 
Feshbach resonance 2l|), the system can exhibit cither 
ferromagnetic 0, ant^ferromagnetic @ , paramagnetic 
24| or even diamagnetic features [25]. 



In this Letter, we investigate the superfluid properties 
of a magnetized spin-I/2 BECs @, \m El]. The mag- 
netization can be obtained (i) from a spin-anisotropic 
condensate prepared or naturally forming in a single spin 
component, or (ii) from a spin-isotropic BEC in the pres- 



ence of an external magnetic field. By showing that the 
spin and density excitations can be completely separated, 
we argue that a spin current, protected against spin ex- 
citations, is possible for a wide range of parameters and 
even in the supersonic regime, where density excitations 
lead to the superfluidity breakdown. This suggests that 
a pure spin superfluid can be regarded independently 
from the usual motional (Landau) superfluid. We test 
the spin-superfluidity criterion by simulating the forma- 
tion and suppression of spin waves in a condensate past 
a magnetic defect and calculate the corresponding drag 
force. 

Spin-l/2 condensate. We consider a spinor BEC with 
two spin projections ±1 (± indices) on the z quantization 
axis. At the mean field level, the system is governed by 
the energy functional 
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in the Madelung representation, where n± and 6± re- 
spectively represent the density and phase of each spin 
component. Here, H = (H x ,H y , H Z ) T is a generic (effec- 
tive) magnetic field. The pseudospin vector S allows the 
mapping to a magnetic system [271 ] and its components 
are linked to the spinor St via the identities 

S x = - (ip+i/j*_+tp* + ip-) = 0hHTcos(A6I) 

S y = l - {W+ - 4'H'+) = N /^Tsin (A0) (3) 



S- 



- (n+ - n_) , 



where A.6 = 6 + — 6- is the relative phase and are 
the intra and intcrcomponent interactions. The S z pro- 
jection can be seen as the magnetization of the system. 
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The free energy of the condensate reads 



F= — fin+(a 1 +a 2 ) — 



- H X S X - H y S y - (H z - AaS z )S z 



(4) 



where n = n + + n_, Aa = a± — a 2 and (i is the 
chemical potential. The last term contains the ap- 
plied magnetic field H z and an intrinsic Zeeman splitting 
i?2 = — AaS z that embodies the spin-spin interaction. 
The ground state can be found from the minimization 
of F, at fixed total density n, with respect to the pseu- 
dospin components and given the normalization condi- 
tion Si + Sy + S 2 = rt 2 /4. Analytic solution can be 
found for the ground state for any H but the expressions 
are cumbersome and it is more instructive to analyze par- 
ticular cases. 

For H z = and H\\ — (H x ,H y ) ^ 0, there are two pos- 
sible solutions. The first one is trivial (S°,S°,S°) T = 
(nH x /2H\\,nH y /2H\\,0) T and corresponds to a param- 
agnetic condensate having its pseudospin aligned with 
the in-plane field Hy = (H x , H y ) T . For a nontrivial so- 
lution with S z 0, which is not the ground state, we 
obtain 
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possible for a spin- anisotropic condensate with Aa =/= 0. 
The pseudospin is aligned with the total effective mag- 
netic field H t = (H x , H y , — AaS z ) T . Interestingly, the 
intrinsic field Hz c = — AaS z u z attributes the conden- 
sate a spin stiffness, preventing the action of Hy. We 
see from Eq.© that S z is characterized by a critical in- 
plane field H c = ±Aan/2 below which the magnetization 
of the condensate would break down towards the trivial 
solution with S z — 0. The chemical potential is found 
minimizing F with respect to the total density n giving 



S z = =$> fj, = (ati + a 2 ) 



H\\, 



(6) 



S z + 

For H z 7^ and H x>y = 0, two solutions are allowed 
as well. The trivial solution corresponds to a magnetized 
condensate with (S°, S°, S°) T = (0,0,±n/2) T . US\\£ 0, 
we obtain 
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As we can see from Eq.([7|, for a given value of H z the 
magnetization S z can change its direction depending on 
the sign of Aa: if Aa > 0, S z is opposed to H = H z u z 
and the condensate behaves as an anti-ferromagnet; in- 
stead, if Aa < the latter corresponds to a ferromagnet. 
Due to the intrinsic field i?zc, the magnetization imposed 
by H z > H c would remain even after switching off the 



applied field. On the other hand, the chemical potential 
reads 

Su = =>■ fji = a±n — H z 
S\\ 7^ =*> fi = (ai + a 2 ) 
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FIG. 1: (Color online) (a) Dispersion of the elementary ex- 
citations of the magnetized condensate. The solid (dashed) 
red line is the positive (negative) Bogoliubov branch of the 
condensed component and the arrows point the formation of 
a Bogolon. The dispersion of the non-condensed component 
(solid blue parabolas) is shifted down for H z < H c or up for 
H z > H c . The wavy line shows the excitation of a plane wave, 
allowed only in the former case, (b) Zero temperature phase 
diagram. SF: Supersonic Flow, SW: Spin Waves formation, 
SS: Spin Superfluid and LS: Landau Superfluid. The dashed, 
solid and dotted blue lines mark the frontier v = c m in the 
cases H z = {— 0.8, 0, 1.2} mc\ respectively. The horizontal 
dashed line stands for v = c s . 

Density and spin waves. Let us now focus espe- 
cially on the case where H x . y = 0, treated as a per- 
turbation in the following, and consider the elementary 
excitations on top of the fully magnetized condensate 
(5°, S°, S°) T = (0, 0, +n /2) T . The minimization of the 
energy functional E leads to the following set of coupled 
Gross-Pitacvskii equations: 

iH ~tt~ = ~ 2~b A ^ + ~~ + ai '^ + | 2 ^+ 

+ a^-l 2 ^-^^, (9) 

dip h 2 2 

ifr—gj- = - —AV>_ - mV'- +"11^-1 V>- 

+ a 2 |V+| 2 V- + : Y^-- (10) 

We start from the Bogoliubov ansatz Q 

il> + = A + e l{ ^ r - ut) + Ble- 1 ^-"*^ 

V>_ = A-e i( - k - r - ut \ (11) 

that involves two counter-propagating plane waves in the 
condensate characterized by the weak amplitudes A + 
and B + , the wavevector k and frequency uj. We seek 
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for a plane wave solution of amplitude A- in the non- 
condensed component. Injecting (fTTj) into Eqs. (|9I1Q[) . 
we obtain the chemical potential /i = ain — H z and the 
following modes 



hu\ L ' 2> = ±y / E(k) (E{k) + 2fi) 
hw„ = E(k) + (2H Z - Aan ), 



(12) 
(13) 



where E (k) — h 2 k 2 /2ni. The elementary excitations on 
top of the condensed component are the well-known den- 
sity waves (or bogolons) following the dispersion uj± = 
±c s k in the long wavelength limit, with c s = \J '/i/m be- 
ing the speed of sound. The static condensate is stable 
against the creation of these excitations, that have a posi- 
tive (| A + 1 2 — \B + \ 2 )hiu^ contribution to the energy, and 
is therefore supcrfluid in the Landau picture. For the 
other spin component, we obtain a parabola shifted by 
a quantity 2H Z — AanQ due to the presence of the con- 
densate (see FigUJa)). The coefficients A + and B + can 
be explicitly found provided the normalization condition 
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where AE = E(k) + fj, - huj ( + } and A^*' = 0. We 
then observe that the excitation of plane waves is for- 
bidden unless the parabolic branch (fTB"]) is redshiftcd 
with respect to the condensate energy. Excitation (by a 
spin-coupling perturbation) of a plane wave Sip- (x, t) = 
A- cxp(ik-r — ut) leads to a spatio-temporal modulation 
of the relative phase A9, leaving the condensate's den- 
sity homogeneous. This kind of excitation can therefore 
be seen as a spin wave, in full analogy with magnetic 
systems (20j , describing the spatial Larmor precession of 
S about the z-axis. The corresponding deviation of the 
pscudospin follows from Eqs. (J3]t4j) , reading 



(^TAE), 



(14) 
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Interestingly, the transverse spin components b~S XtV oscil- 
late out of phase, while the magnetization SS Z remains 
unchanged. This simply means that spin and density 
waves are decoupled. Moreover, the parabolic branch 
(fTB"]) is blue-shifted in respect to the anti-ferromagnetic 
condensate if H z > H c = Aan^/2. The condensate is, in 
that case, stable against the creation of both density and 
spin waves, guaranteeing the simultaneous occurrence of 
Landau and spin superfluidity. We note that the ferro- 
magnetic condensate automatically fulfills this require- 
ment, even if H z = 0. 



To fully understand the physical meaning of spin su- 
perfluidity considered in this Letter, we should address 
the situation of condensate propagating with a velocity 
v. In that case, the dispersions (|12j|13p arc shifted by 
the quantity Hv ■ k. According to the Landau's criterion, 
the superfluidity is lost for v > c s = y^x/m when the 
creation of bogolons becomes energetically favorable. In 
the same spirit, the formation of spin waves becomes pos- 
sible as soon as the branch of spin waves acquires some 
negative energy (see insets of Figj3]), allowing the energy- 
conserving transfer from the ui^P to u;_ branches. Such 
process occurs above the critical velocity 



AH Z - 2Ac 




(16) 



The above expression allows us to construct the phase 
diagram shown in Figfljb), containing four regimes (see 
captions): exclusive Landau (c m < v < c s , pink region) 
or spin (c s < v < c m , purple region) superfluidity, the 
absolute superfluidity (v < c s ,c m , light blue region) and 
the onset of both spin and density waves (v > c s ,c m , 
white region). 

Generation of spin waves. The most famous exper- 
iment to highlight superfluidity involves a condensate 
propagating against a potential barrier. As it is known, 
in the supersonic regime v > c s , backscattered density 
waves appear upstream from the defect; in the subsonic 
regime, however, no radiation occurs. 




(15) FIG. 2: (Color online) (a) Spin waves excited by a magnetic 



defect. The S x , S y and S z components are shown with red, 
blue and black lines respectively. The dashed, dotted and 
solid lines correspond to U m /Hc = {0.1, 0.4, 0.7} respectively, 
(b) Spin drag force Fd experienced by the condensate. The 
dashed black, dashed dotted blue and solid red line corre- 
sponds to H z = {0,1.4, 0.5} mc 2 8 and a 2 = {1.2,0,1.2} re- 
spectively. 

It is therefore natural to investigate what happens to 
the flow if a magnetic defect H x (r) is placed instead of 
a usual potential. For definitcness, we take H z = and 
Aa > in the following discussion. A first expected 
effect regards the conversion from the + to the — com- 
ponent via spin rotation. Naturally, the conversion effi- 
ciency depends on both the strength of the field and on 
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the speed of the flow. Second, beyond the defect, the 
residual magnetization Hz c leads to self-induced Larmor 
precession of the pseudospin, thus creating a spin wave as 
discussed previously. To get some intuition we can con- 
sider a zero-dimensional system in which the magnetic 
defect H x is switched on and off. In such configuration, 
the dynamics of the pseudospin can then be described by 
an undamped Landau-Lifshitz equation hdtS = H x S, 
which leads to the following set of nonlinear equations 

hd t S x = AaS z S v , (17) 
hdtSy = (-AaS x - H x ) S z , and (18) 
hd t S z = H x S y , (19) 

to be solved with respect to the initial condition So = 
(0,0, no/2). Here, we assume that H x = H x (t)u x with 
H x {t) = U m exp[— (t — t ) 2 /At 2 ] representing a pulse of 
duration At switched on at the instant to- Eqs. (|17I 
[T9l) can be solved analytically in the subcritical regime 
U m <C Hzc- Up to a small phase factor, the solution 
reads 

S x (t) = --£Ls°em(H$j) (20) 
S y (t) = +^S° y cos (Hit) (21) 

s -«> - \/?-(^) 2 < 22 > 

at long times, where H% c = AaS z , At m = 
Hn/y/U* + H'l c and 5° z are given by Eq.©. Full nu- 
merical solutions to Eqs. fj 1 741 1 9 [) are presented in the 
FiglJJa), showing the dynamics of the pseudospin pre- 
cession. We can observe the formation of out-of-phase 
oscillations of S x and S y . By increasing the value of U m , 
S z (and, consequently, 7?ze) is reduced, which in turn in- 
creases the magnitude and reduces the frequency of the 
spin wave as captured by Eqs. (|20H2"2")) . 

To further understand the differences between the 
usual and spin superfluidity, we address the question 
of dissipation in the spin flow. In the presence of a 
mangnetic defect, the drag force is defined as F m = 
— / ('0+ V-fTr'0- +c.c)dr. For simplicity, we consider a lo- 
calized defect of the type H x = U m 8(r), such that F m = 
Ud J"(V>+VV>-|r=o+ce)eJr. Expanding the wave functions 
according to (TTTj) , and using the Hopfield amplitudes (fT4")) 
modified by the Doppler shift ftk-v, we obtain, in first or- 
der in A± and B + , F m = U rn n j ikSf(k, k ■ v)dk/(27r) 2 , 
where 

i p 

S f (k,uj) = —- — r (23) 

JV ; 2m(co + iQ) 2 - uPl y ' 

represents the dynamic structure factor associated with 
the spin excitations. The infinitesimal positive imaginary 
part +i0 is added to the frequency lu to be consistent with 



the usual Landau causality prescription. Computing the 
integral of ((23)) in the complex plane by means of the 
relation 1/z + iO = P(l/ z) + iirS(z), and noticing that the 
integral over ip (<p is the angle between k and v) vanishes 
for the principal part P as a consequence of symmetric 
integration limits, we easily obtain the explicit expression 

F »=^(" + V^)V (24) 

We observe that the force is zero for v < c m and mono- 
tonically increases with v if v > c m (see FigJ^b)). For 
overcritical flows, v 3> c m , the force reduces to the 
classical result F m cx v. Notice that Eq. (|24[) is for- 
mally similar, but physically different, to the drag force 
Fd = a 2 d n m 2 [v 2 — c 2 )/(H 3 v)u v acting on the superfluid 
density in the presence of a localized defect Vd = adS(r) 
[28| . where the role of critical velocity c m is played by c s . 

Numerical results. To illustrate the results discussed 
previously, we performed numerical solutions of Eqs. 
(|9ll0p for parameters corresponding to different regions 
of phase diagram depicted in FigOJb). We considered 
a 2D homogeneous magnetized condensate propagating 
with the velocity v = v x u x , by setting the condition 
*o ( r ) = {\f n He- 1 ' kx 1 0) T , with k x = m/hv x . In order 
to distinguish between the Landau and the spin super- 
fluidity, we plugged a potential barrier and a magnetic 
defect, well separated from each other, to respectively 
excite density and spin waves in the flow. The latter 
demands an extra local coupling term —H x (r)ip^z/2 in 
the equations. The results are summarized in Fig{3j 
where the four regions of the phase diagram are consid- 
ered. We enclose snapshots of the linear wave function 
i>x = tp+ + ip- j which is the best suited quantity to il- 
lustrate density and spin waves simultaneously. In each 
panel, the insets show the corresponding dispersion of 
elementary excitations. Panel (a) corresponds to a su- 
personic flow configuration in which spin waves can de- 
velop as well (u > c s , c m ). Interestingly, while the density 
waves form upstream from the barrier, the spin waves are 
found downstream due to their positive group velocity 
Panel (b) shows the exclusive spin superfluid regime in a 
supersonic flow, where only density excitations take place 
(c s < v < c m ). Panel (c) illustrates the exclusive Landau 
superfluidity, for which only spin waves can be excited 
(c m < v < c s ). Finally, panel (d) demonstrates the abso- 
lute superfluidity of the spinor condensate, in which no 
excitations of any kind are observed (v < c s , c m ). 

Conclusions. We have shown that magnetized two- 
component condensates allow the formation of both den- 
sity and spin waves that can be suppressed independently 
depending on the type of spin-spin interaction, the ap- 
plied magnetic field, and the speed of the flow. In such 
system, superfluidity requires both the Landau and spin 
superfluidity criteria to be fulfilled. We have constructed 
the phase diagram of the system and illustrated the lat- 
ter with a propagating condensate against a regular and 
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FIG. 3: (Color online) The four phase diagram regimes illus- 
trated numerically with both a regular (upper) and a mag- 
netic (lower) defect within the same flow. The colormap shows 
\ipx\ 2 and the insets display the dispersions. The visual order 
of the phases is preserved with respect to Fig[ljb). 



magnetic defect and calculated the spin drag force. A 
well suited system to implement experimentally our pro- 
posals would be a circularly polarized exciton-polariton 
condensate either in ferromagnetic regime [2j| or under 
an applied magnetic field to compensate for the antifcr- 
romagnctism [251 130( . This system is to our knowledge 
the one offering to tune the parameters of a flowing con- 
densate in the most flexible way [3l[. Finally, the obser- 
vation of superfluid propagation of spin currents would 
represent an important milestone in spintronics and spin- 
optronics. 
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